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ENERGY SPECTRA OF ELEMENTS WITH 18 < Z < 28 
BETWEEN 10 AND 300 GeV/amu 

Michael D. Jones*, J. Klarmann*, E. C. Stone^, C. J. Waddington c , W. R. Binns a , 

T. L. Garrard , M. H. Israel* 

a) Washington University, St. Louis, MO 63130, USA 

b) California Institute of Technology, Pasadena, CA 91125 USA 

c) University of Minnesota, Minneapolis, MN 55455, USA 

1. Introduction . The HEAO-3 Heavy Nuclei Experiment (Binns, et al., 1981) is 
composed of ionization chambers above and below a plastic Cherenkov counter. We 
have measured the energy dependence of the abundances of elements with atomic 
number, Z, between 18 and 28 at very high energies where they are rare and thus need 
the large area x time of this experiment. We extend the measurements of the Danish- 
French HEAO-3 experiment (Englemann, tt a L, 1983) to higher energies, using the 
relativistic rise of ionisation signal as a measure of energy, and determine source 
abundances for Ar and Ca. 

2. Data Analysis . We confine this analysis to events in the 1.1 m 2 sr in which the 
cosmic rays penetrate all six ionization chambers, and to the first 314 days of the flight, 
when all six ionization chambers were functioning properly. These selections give the 
highest possible ionization resolution. We select only events incident with geomagnetic 
cutoff greater than 8 GV, so Z is determined by the Cherenkov signal. We further 
require agreement between the means of the three ionization chamber signals above the 
Cherenkov and those below, to eliminate most nuclear interactions inside the 
instrument. 



Figure 1 



Figure 2 


Figure 1 is a schematic plot showing the locus of events for each element, 20 < Z 
< 28. Zc is the square-root of the Cherenkov signal normalized so Zc ~ Z at high 
energy. Zi is the square-root of the ionization normalized so Zi = Z at minimum- 
ionizing. Figure 2 is a histogram of Z: for events with 29.3 < Zi < 29.5, the region 
between the dashed lines in figure 1. This histogram includes Fe, Co, and Ni at about 
130, 34, and 12 G< V/amu respectively. The abundance of each element in each of 
eighty such histograms is determined by maximum-likelihood fitting. 
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Corrections to these raw abundances were calculated to account for interactions 
both in the lid in front of the first ionization chamber and in the Cherenkov counter 
and other material between the ionization chambers. The latter calculation included 
only charge changes of one or two charge units, other interactions having been 
eliminated by the requirement for agreement between upper and lower ionization 
chambers. The calculation assumed that at all energies the total cross-sections were 
given by the formula of Westfall et aL (1979) and the partial cross-sections for changing 
by n charge units were the same fraction of the total cross-section for any projectile as 
was measured by Webber and Brautigam (1982) at 980 MeV/amu for Fe on C. These 
calculated corrections lowered the raw abundance ratios by typically 10 to 30 percent. 

3. Energy Scale. We used Zi/Z as a measure of energy, and derived an empirical 
energy calibration by comparing our Fe observations with an Fe energy spectrum 
derived from a compilation of previously published measurements (Webber, 1982). This 
Fe spectrum was multiplied by an empirical geomagnetic transmission function which 
represented the fraction of time when the geomagnetic cutoff permitted Fe of that 
energy to reach the instrument; the product was the effective Fe spectrum at the 
instrument, averaged over many orbits. This energy spectrum was then converted to a 
Zi/Z spectrum using a trial form of the energy dependence of Zi/Z. Finally this 
calculated spectrum was folded with the instrument’s ionization resolution, and the 
resulting Zi/Z spectrum was compared with the data. The process was iterated, by 
changing the assumed form of the energy dependence of Zi/Z, until the calculated and 
observed spectra of Fe agreed. 

The resulting energy dependence is consistent, between about 10 and 100 
GeV/amu, with one derived independently for a different detector system by Barthelmy 
et cl. (1985, OG 4.1-7). Above about 200 GeV/amu, the shape of our calibration curve 
depends upon the assumption we made that the Fe energy spectrum falls as a power- 
law with exponent -2.7 at energies above those where it has been measured. 

Abundance ratios derived from data in a particular Zi/Z histogram were plotted 
at the mean energy from which those particles came, assuming the Fe energy spectrum 
and the calibration curve described above. The energy resolution implied by our 
ionization resolution of 0.40 charge units is comparable to the spacing of the points in 
figure 3. 

4. Results. Figure 3 gives the resulting abundances of several elements relative to Fe 
as a function of energy. The X symbols are the results of the DF experiment 
(Englemann, 1983), while the O symbols are the results of this experiment. In every case 
o-ir resuits are consistent with those of the DF experiment in the interval where both 
experiments apply, 10 to 25 GeV/amu. At higher energies our data generally continue 
the DF trend. 

Our data above 10 GeV/amu suggest a Ni/Fe ratio slightly dependent upon 
energy, with a best fit power law of exponent -0.050 ± 0.016. If we igno-e this slight 
variation with energy, then the mean vaiue of the Ni/Fe ratio over our data is 0.054 ± 
0 . 001 . 

For the secondary ratios, K/Fe, Sc/Fe, Ti/Fe, and V/Fe, our data indicate an 
extension to about 100 GeV/amu of the same power law dependence a; was indicated 
*\v the DF data. Figure 4 shows the best fit exponent for each of these ratios combining 
the DF and our data. The steepening of the slope with decreasing Z is expected as 
lower Z elements have greater contributions from tertiary nuclei. 








9J 01 3AI1V13H 30NVQNnaV 


* 


100 







4 


OG 4.1-8 


- 1 - 


Our data indicate a leveling of the Ca/Fe ratio above the energies of the DF 
experiment, as would be expected from an energy-independent primary component 
becoming increasingly significant at higher energies as the secondary component 
becomes less abundant. We fitted the combination of the DF and our ratios to a 
function aE p +t>. With p = -0.2 J5 ± 0.010, interpolated from figure 4, the result is 
primary Ca/Fe -* 0.094 ± 0.004 A galactic propagation calculation on the DF data 
(Lund, 1984) gives a source abundance of Ca/Fe = 0.065 ± 0.019. 



Figure 4 

Exponent of power law fit to 
abundance of element Z 
relative to Fe. 


A similar fit to the energy dependence of Ar/Fe, but with p = -0.33 ± 0.01, gives 
primary Ar/Fe = 0.026 ± 0.003. Propagation of the DF data (Lund, 1984) gives Ar/Fc 
= 0.032 ± 0.008. Source abundances inferred from such galactic propagations on 
observed abundances 2.5 to 5 times higher must depend critically upon the adopted 
fragmentation cross-sections, while our extension of observations above 100 GeV/amu 
permits inference of source abundances without galactic propagation calculations. 

5. Acknowledgement. This work was supported in part by NASA grants NAG 8-498, 
500, .502. and NCR 05-002-160, 24-005-050, and 26-008-001. 
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Elemental Abundances of Cosmic Rays with Z > 33 as Measured 

on HEAO-3 

BJ. Newport 0 , E C. Store* 1 . C.J Waddington b , W.R. Binns c . T L. Garrard J . 

( M.H. Israef, and J. Klarmanrf 

California Institute of Technology. Pasadena, California 91125, USA 
University of Minnesota. Minneapolis. Minnesota c 5455, USA 
c Wash;ngton University, St Louis. Missouri 63130, USA 


1. Introduction 

The Heavy Nuclei Experiment on HEAO-3 (Binns et al., 1981) uses a combination 
of ion chambers and a Cerenkov counter. During analysis, each particle is assigned two 
parameters, Zc and Z[, proportional to the square roots of the Cerenkov and mean ioniza- 
tion signals respectively. Because the ionization signal is double valued, a unique assign- 
ment of particle charge, Z, is not possible in general. Our previous work (Binns et al., 
1983, 1985, and Stone et al., 1983) has been limited to particles of either high rigidity or 
low energy, for which a unique charge assignment was possible, although those subsets 
contain less than 50% of the total number of particles observed. In this paper we discuss 
the use of the maximum likelihood technique to determine abundances for the complete 
data set from —1.5 to —80 GeV/amu. 

Figure 1 shows the possible values of Zc and Z\ for elements rear iron, and indicates 
the substantial overlap between adjacent elements, even before smearing by the resolution 
function. In Figure 2, the curves of Figure 1 have been transformed using the variable 
Zc/Z| instead of Z\. This transformation simplifies the following data analysis. 

2. Analysis 

Panicles were selected from the full exposure. 580 days, and were required to have a 
good Cerenkov signal, at least one good ion chamber and a reliable trajectory. These par- 
ticles were assigned an initial charge estimate. Z**. and 1/40 of those wiih Z^ >19.5 were 
saved, together with all the remaining particles with Z„,>30. The selected panicles were 
binned in a two dimensional histogram, with one axis being the logarithm of Zc and the 
other being Zc/Z|. Figure 3 shows a contour plot of the region of this histogram near 
iron. 
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We have used the iron distribution as a reference distribution for the other elements, 
by scaling it according to the following criteria: 

1) All elements are assumed to have the same energy spectrum, although it is known 
that the sub-iron secondary elements have steeper spectra (e.g. Jones et ai., 1985, OG 
4.1-8). 

2) Energy independent scaling factors have been used. 

3) The resolution of the instrument is a constant fraction of the signal. 

4) Non-Z 2 corrections to the scaling laws for Zc and Zc/Z| have been determined 
directly from the data. 

Those elements contaminating the iron distribution have been approximately 
removed by scaling the contaminated distribution and subtracting according to an 
assumed set of abundances. The resulting "clean" iron distribution was smoothed and 
then scaled to the high Z elements, using cubic interpolation techniques. 

The likelihood of a given set of abundances may be calculated using Poisson statis- 
tics, and maximized by iterating until all the first derivatives are zero, yielding the best fit. 
Results from the application of this method will be reported. 

3. Acknowledgements 

We thank B.W. Gauld for assistance in programming for data analysis. This work 
was supported in part by NASA grants NAG 8-498. 500, 502 and NGR 05-002-160. 24- 
005-050, 26-008-001. 
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ABUNDANCES OF ‘SECONDARY’ ELEMENTS AMONG THE 
ULTRA HEAVY COSMIC RAYS - RESULTS FROM HE AO- 3 

J. K!armann\ S. H. Margolis 1 , E. C. S' one\ C. J .Waddington', 

VV. R. Binns*, T. L. Garrard 11 , M. H. i rael* and M. P. Kertzman' 


‘Department of Physics and the McDonnell Center for the Space Sciences, 
Washington University, St. Louis MO 63130, USA. 
b George W. Downs Laboratory, California Institute of Technology, 

Pasadena CA 91125, USA. 

c School of Physics and Astronomy, University of Minnesota, 

Minneapolis MN 55455, USA. 

1. Introduction. This paper discusses observations of the abundances of elements of 
charge 62<Z<73 in the cosmic radiation from the HEAO-3 Heavy Nuclei Experiment 
(HNE). These elements, having solar, and presumably sourct, abundances much less 
than the heavier Pt and Pb groups, are expected to be largely products of spallation. 
Thus they are indicators of the conditions prevailing during the propagation of cosmic 
rays. The abundances have changed from those reported previously (Klarmann et al., 
1983) due to a different data selection (Binns et al., 1985). This resulted in better charge 
resolution and in a higher mean energy for the particles. All the particles we have 
included in this paper were required to have had a cutoff rigidity R c > 5 GV. This 
allowed the charge determination to be based solely on the Cherenkov measurement. 
For a description of the detector see Binns et al., (1981). 

3. Analysis. The data selection in this paper is identical to that of Waddington et 
al., (1985, OG4.4-7). We have considered only the following physically significant groups 
of charges: 


Name 

Abbreviation 

Range 

Number observed 

Lead and Platinum 

PbPt 

74< Z<86 

52 

Heavy secondary 

HS 

70< Z<73 

10 

Light secondary 

LS 

62<Z<69 

34 


Our discussion will be in terms of the ratios: HS/'PbPt and LS/PbPt. In the table, 
column a) shows the results observed in the detector. The correction factor to outside 
the detector was derived by propagating eight different plausible theoretical abundances 
outside the detector through slabs of hydrogen approximating the distribution of alumi- 
num traversed by the particles going into and through the detector. The change of the 
abundance ratios from outside the detector to inside was nearly independent of the origi- 
nal ratios and is given as a multiplicative coirection factor in column b). The abundance 
outside the detector, column c) is the product of columns a) and b). 


Ratio 


IEAO Results 

Ariel 

HEAO/Ariel 

Inside 

Detector 

a) 

Correction 

Factor 

b) 

Outsid; 

Detector 

c ) 

Outside 

Detector 

d) 

O >(de 
Dt .ector 
e ) 

HS/PbPt 

LS/PbPt 

0.19±0.07 

0.65±0.!1 

0.85±0.02 
0.87 ±0.02 

0.16±0.06 
0.57±0. 12 

0.27 ±0.07 
0.88±0. 15 

0.59±0.27 
0.65±0. 18 


r. 


Results from the Ariel-0 UH-nuclei detector which was exposed in a 55 * inclination 
orbit (Fowler et al., 1984) are given in column d), while column e) gives the ratio of our 
HEAO results to those of Ariel. It is seen that for both ratios our result is about 60 % 
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to 65% that of Ariel's. While these differences are only significant at a level of 1.5 to 2.0 
standard deviations, it is unlikely that they are just statistical fluctuations. The data of 
Ariel extend to significantly lower energy than ours. At lower energies the abundance of 
secondaries is expected to be greater since both the interaction cross sections and the 
escape length are latger. We cannot tell yet whether this energy dependence is sufficient 
to explain the difference. 

3. Comparison with Models. The abundance ratios can be compared with predic- 
tions of various models. The source abundance used was either the solar system abun- 
dances of Anders and Ebihara (1982) (No FIP) or the same adjusted for an exponential 
dependence (Brewster et al., 1983a) on the first ionization potential (FIP). These were 
then propagated through the interstellar medium, assuming a leaky-box model, and 
using the revised code of Brewster et al., (1983a, 1985) with a rigidity dependent escape 
length (Ormes and Protheroe, 1983) that is 6.21 g/cm of hydrogen at 7 GV. The calcu- 
lated values are for approximately the same mix of rigidities as the HEAO data. A 
different model of FIP fractionation (Cook et al., 1979; J. P. Meyer, 1981), in which the 
cosmic ray source is suppressed by a constant factor relative to solar abundances for ele- 
ments with ionization potential above 9 eV, yields propagated abundance ratios which 
in, this charge range, are indistinguishable from those of the unfractionated source. 
Similarly, propagation of an r-process source abundance yielded ratios which in this 
charge region were close to those from a solar system source. Neither of the last two 
results is plotted in figure 1. 



Figure 1: Comparison of the observed and predicted abundance ratios. 

In the No FIP alternate’ propagation an independent code was used (Margolis, 1983) to 
predict abundance ratios after propagation through leaky boxes of various escape 
lengths. The results were then combined using the same rigidity dependent escape 
length distribution as above to yield the inverted triangle point in figure 1. With this 
rigidity de^rendent distribution the mean escape length encountered by the observed par- 
ticles is ~3g/cm : . This point, when compared to the other No FIP point, is an 
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indication of the variation possible in the propagation calculation. The point labeled 6 
- /cm-*’ in figure 1 is the result of the same propagation through a leaky box with a sin- 
gle escape length of 6 g/cm 2 of hydrogen. The difference between this point anc the 
‘alternate’ point shows the dependence of the results on the escape-length distribution. 
In figure 1 experimental values are solid with error bars. 

The dependence of the abundance ratio on propagation can also be demonstrated in 
a different way. Every point in figures 2 and 3 (Margolis and Blake, 1985) corresponds 
to the calculated ratio after propagation of a solar system source without FIP through 
hydrogen with a mean free path distribution rising linearly from zero to the desired 
‘truncation’ then falling exponentially with the given ‘escape length (Margolis, 1983). 



ESCAPE LENGTH 


Figure 2: LS/PbPt 



ESCAPE LENGTH 


Figure 3: HS/PbPt 

In contrast to the data in figure 1, these figures assume that all particles tra-.ers; the 
same path length distribution. Our results are represented by the cross-hatched region 
yielding possible combinations of escape length and truncation. 
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A a expected the predictions in this charge region are nearly independent of escape- 
length since the interaction mean free path is so short. However the results do not agree 
with more than a minute amount of truncation of short path lengths. The fact that at 
zero truncation an escape length of ~2g/cm J is indicated seems to support the rigidity 
dependent escape length proposed by Ormes and Protheroe, (1983). 

4. Discussion. Our observed values of the secondary ratios are in reasonable agree- 
ment with the prediction based on a model without FIP fractionation or with a step 
function FIP fractionation at the source; however, our observations are in distinct 
disagreement with the models that include exponential FIP fractionation. This is con- 
trary to the conclusions found at lower charges (Binns et a!., 1982, 1983) where observed 
abundances agreed better with those expected from a solar system source with FIP frac- 
tionation than without. Thus other representations of source fractionation may be 
involved. 

Our results do fit the predictions obtained using the standard leaky box model in 
this energy range. The applicability cf this model to lower energies requires further 
investigation. 

5. Acknowledgement*: This work was supported in part by NASA giants NAG 

8-448, NAG 8-498, 500, 502 and NGR 05-002-160, 24-005-050, 26-008-001. 
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LEAD, PLATINUM, AND OTHER HEAVY ELEMENTS 
IN THE PRIMARY COSMIC RADIATION— HEAO-3 RESULTS 

C.J. Waddington, 1 W.R. Binns, 2 N.R. Brewster, 1 D.J. Fixsen, 1 
T.L. Garrard, 3 M.H. Israel, 2 J. Klarmann, 2 B.J. Newport, 3 E.C. Stone 3 

School of Physics and Astronomy, University of Minnesota, 

Minneapolis, MN 55455 

department of Physics and the McDonnell Center for the Space 
Sciences, Washington University, St. Louis, MO 63130 
3 George W. Downs Laboratory, California Institute of Technology, 

Pasadena, CA 91125 

1. Introduction. This paper reports an observation of the abundances of 
cosmic-ray lead and platinum-group nuclei using data from the HEAO-3 
Heavy Nuclei Experiment (HNE) which consisted of ion chambers mounted on 
both sides of a plastic Cherenkov counter (dinns et al., 1981). 
Previously we have reported on a search for actinide nuclei, Z > 88 
(Binns, et al. 1982a). Further analysis with more stringent selections, 
inclusion of additional data, and a calibration at the LBL Bevalac, have 
allowed us to obtain the abundance ratio of lead and the platinum group 
of elements for particles than had a cutoff rigidity R c > 5 GV. 

2. Analysis. We have analyzed 580 days of exposure and considered 
selected data for those events where the Cherenkov detector and at least 
two of the ion chambers were triggered. These selection criteria will 
be described elsewhere, Binns et al. (1985). 

Two sets of events satisfying the selections were formed — one for 
which Z > 49.5; the other, a "normalization" set, with 1/400 of all 

events with Z > 19.5, chosen at random. 

The events were separated into two groups, 67Z with R^ > 7 GV and 
33Z with 5 < R <7 GV. The charge scale and resolution for each group 
were determined independently by examining the iron peak in the 
corresponding normalization set. 

In both groups, the nuclear 
charge of each event was 50 
inferred from the Cherenkov 
signal, assuming that the signal 
was simply proportional to Z 1 , £ 

Garrard et al. (1983). § 20 

Fig. 1 shows the observed * 
charge spectrum. This data set £ 
demonstrates an odd-even £ 
abundance effect for 5QsZs56 and § „ 

.1 sharp falloff in abundances 
between 5b and 60, similar to 
that found previously in a data 
subset having higher charge 
resolution (Binns et al. 1983). 0 

The 322 nuclei with Z 2 50 used 
in this p alysis correspond to 

( 9.6 t 0.5)10* iron nuclei which Fig. 1. Observed charge spectrum 
satisfy the same selection with charges assigned assuming a Z 2 
criteria and are observed within dependence of the Cherenkov signal. 
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the instrument, not in free space. The quoted uncertainty is 
predominately due to the uncertainty in resolving 25 ^ from 26 ^ e> 

3. Comparison with Other Data. Results that cover this charge range 
have been reported from the Ariel -6 UH-nuclei detector which was exposed 
in a 55° inclination orbit (Fowler et al. 1984), and hence extends to 
appreciably lower energies than our data. In order to analyze the Ariel 
data, Fowler et al. had to deconvolve their data using an extrapolation 
of the resolution function found for Fe and lighter nuclei. We have not 
attempted a deconvolution of our charge spectrum, since the results of 
such a process are quite sensitive to the form of the assumed resolution 
function, particularly when individual element peaks are not apparent in 
the data. Due to our limited charge resolution we have considered only 
the following physically significant groups of charges: 

Name Abbreviation Range Number observed 

"Lead" Pb 81SZS86 10 

"Platinum" Pt 74SZS80 42 


The ratio of the abundance of lead to platinum will be compared 
with other data and with model predictions. The secondary ratios will 
be discussed elsewhere, see Klarmann et al. (1985; 0G 4.4-6). 

The value of 0.24±0.08 for the Pb/Pt ratio derived from our 
observations differs from that outside the detector because of nuclear 
interactions during entry and penetration of the detector and the 
instrumental resolution, which smears the charge distribution. For each 
of eight plausible models we calculated abundances expected near earth, 
as described below. Entry into the detector was then simulated by 
propagation through various slabs of hydrogen approximating the amount 
of aluminum in the various paths into aid through the detector. The 
resulting element distribution inside the detector was then convolved 
with the instrument resolution to derive the distribution we would 
expect to observe. Although the eight models gave very different values 
for the ratio at the outside of the instrument, the factor by which the 
ratio changed after propagation into the instrument and convolution with 
the resolution was nearly the same for all the models. Therefore, we 
have used a single correction factor of 1.06±0.02 for the ratio. 

Our resulting ratio of 0.25±0.09, outside the detector, can be 
compared with the corresponding result reported by the Ariel experiment 
of 0.40±0.10. If this result is combined with those on the secondary 
ratios, there does seem tw oe a systematic difference between the two 
sets of results, although of marginal statistical significance on any 
individual ratio. 


4. Comparison with Models. Our observed charge spectrum, Fig. 1, can be 
compared with those predicted by various models. A series of 
predictions were made using the solar system abundances of Anders and 
Lbihara (1982) and the derived s- and r-process contributions to these 
abundances. These abundances, taken as calculated, or adjusted for the 
effects of an exponential dependence on first ionization potential (FIP) 
fractionation, were used as source abundances. An alternative 
dependence on FIP, with l step at 9 eV or above, Cook et al. (1979); 
Meyer (1981), would lead to abundances essentially independent of FIP. 
These various source abundances were then propagated through the 
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interstellar medium, assuming a leaky-box model, and using the revised 
code of Brewster et al. (1983 , 1985) with a rigidity dependent escape 
length (Ormes and Protheroe 1983) that is 6.21 g/cm 1 of hydrogen at 7 
GV. We have used the cross-sections calculated from the formalism of 
Silberberg and Tsao (1983). The predictions of this program are in good 
agreement with the latest predictions obtained by Margolis and Blake 
(1983), at least for the solar system source abundances. 

In Fig. 2, we have shown the calculated values of this ratio for 
solar system abundances and for r-process abundances; s-process 
abundances are not given because they show little relation to the 
observed values with Pb/Pt ratios of n, 1.0. 


Fig. 2. The "lead to 
platinum" ratio as observed 
and predicted. Observed 
values are shown shaded, 
while in space values are 
shown solid and with error 
bars. The shaded and solid 
triangles indicate the 
ratios when a non Z 1 
correction to our charge 
assignments is included. 

Our observed ratio for Pb/Pt (Fig. 2) is distinctly lower than thct 
predicted from solar system source abundances in any of the four models 
considered. In particular, even considering the models without 
exponential FIP fractionation, we find an observed ratio that is 
distinctly lower than that predicted for either a solar-system or an r- 
process source. This result might suggest that, unlike the cosmic rays 
with Z < 60 (Binns et al. 1982b, 1983), the cosmic rays with Z around 80 
come from a source with a distinctly different nucleosynthesis history 
Lhan do the solar system elements. However, two alternatives to this 
conclusion must also be considered. First, the Pb abundance in the 
cosmic ray source may be suppressed by some form of source fractionation 
which depends upon a different parameter chan FIP. Second, it could be 
that the Pb abundances assumed in our model calculations are not really 
representative of the solar system or of the r- or s- process 
contributions to the solar system. 

We have noted (Israel et al. 1983) that the cosmic ray abundance of 
Ge relative to Fe is down by a factor of about two compared to the solar 
system. Ge, like Pb, is one of the fpw volatile elements with moderate 
to low FIP. The factor-of-two underabundance of Ge lends support to the 
suggestion (Cesarsky and Bibring 1980; Epstein 1980) that it is volatile 
elements, rather than elements with high FIP, which are underabundant in 
the cosmic rays. Such a source fractionation dependent on volatility 
could produce our observed low Pb abundance even with a cosmic ray 
source whose composition is essentially the same as that of the solar 
system. 

Alternatively, there are reasons for believing that the source 
abundances of Pb used in our models may not be representative of the 
solar system values. Our observed Pb/Pt ratio could be consistent with 


I I >--->■> » I > I I 

r- proctii 
SH«r 

MIC1. 

HC AO 

I I » I * * » I ■ « — ♦— * — 

o o .2 a« a« aa 

LEAP / PLATINUM 



* 

I 



4 


OG 4.4-7 


that expected from a "Pb-poor r-process", either with or without FIP 
fractionation. 

It is possible that the assumed solar system Pb abundance itself is 
too high. If the Anders and Ebihara Pb abundance were twice that of 
typical solar system matter, then a solar system source abundance, 
either with or without FIP fractionation, would agree with our data. 

Finally, we note thac Ge and Pb, like most elements with higher 
FIP, have abundances in C2 chondritic meteorites about a factor of two 
lower than abundances in the Cl chondrites which are the basis for the 
Anders and Ebihara solar system abundances. If the C2 rather than the 
Cl chondrites were more nearly representative of the composition of the 
heavier elements in the solar system, then our low Pb/Pr ratio would 
again be consistent with a cosmic ray source of composition similar to 
that of the solar system. 

Thus, while our Pb/Pt ratio is distinctly lower than that predicted 
by any of the standard models for cosmic ray sources, it is possible 
that the difference is not an indication that the cosmic ray source 
composition is greatly different from that of the solar system, but 
rather that there is less Pb in the solar system and in the r-process 
than is assumed in the standard model. 

5. Acknowledgements . The research was supported in part by NASA under 
grants NAG 6-498, 500, 502, and NGR 05-002-160, 24-005-050, and 26-uu8- 
001 . 

References 

Anders, E. , and Ebihara, M. (1982), Geochimoca et Cosmochimica Acta , 46, 
2363. 

Binns, W. R., et al. (1981), Nucl. Inst. Meth. , 185, 415. 

Binns, W. R., et al. (1982a), Ap.J. (Letters) . 261, LI 1 7 . 

. (1982b), Ap.J. ( Letters ). 24 7, LI 15. 

Binns, W. R., et al. (1983), Ap.J. (Letters ). 267, L93. 

Binns, W. R., et al. (1984), Adv. Space Res. . 4, 25. 

Binns, W. R., et al. (1985), to be published in Ap.J. , Oct. 1. 

Brewster, N. R., Freier, P. S., and Waddington, C. J. (1983), Ap.J. , 
264, 324. 

. (1985), to be published in Ap. J. , July 1. 

Cesarsky, C. J., and Bibring, J. P. (1980), IAU Symposium , 94, 361. 

Cook, W. R., et al. (1979), 1 6th ICRC . Vol. 12, 265. 

Epstein, R. I. (1980), MNRAS . 193, 723. 

Fowler, P. H. , et al. (1984), 9th European Cosmic Ray Symposium . 

Garrard, T. L. , et al. (1983), 18th ICRC . Vol. 9, 367. 

Israel, M. H., et al. (1983), 18th ICRC . Vol. 9, 305. 

Klarmann, J., et al. (1985), this conf., OG 4.4-6. 

Hargolis, S. H., and Blake, J. B. (1983), 18th ICRC, Vol 9, 283. 

-leyer, J. P. (1981), 17th ICRC . Vol. 2, 281. 

Silberberg, R., and Tsao, C. H. (1983), Ap.J. Supp. 25, 315, 335. 




r. 


l 


^ 'A 


N85-32049 

1 OG 7.1-4 

IMPLICATIONS Cc SOURCE ABUNDANCES OF ULTRAHEAVY COSMIC RAYS 
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S. H. Margolis 4 , E. C. Stone*, C. J. Waddington e , 

4 Washington University, St. Louis, MO, USA 
* California Institute of Technology, Pasadena, CA, USA 
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1. Introduction . In this paper we will examine the ratio of cosmic-ray source 
abundance to solar-system abundance for individual elements. In particular we will 
look at correlations of these ratios with first-ionisation potential (FIP) and also with 
the expected mass-to-charge ratio (A/Q) of the elements in a million-degree plasma. 
We have previously examined the FIP correlation and shown that the correlation is 
affected by the choice of C2 or Cl chondritic meteorites as the solar-system standard 
for comparison (Binns, et a/., 1084). An A/Q correlation was suggested by Eichler and 
Hainebach (1081), as a consequence of their model of shock acceleration in the hot 
interstellar medium, and has been examined by Israel (1085). These correlations are 
presented in the following four figures. 

2. Explanation of F**ures . Figure 1 plots the ratio of cosmic-ray source to solar- 
system abundances, normalised to unity for Fe, as a function of the first ionisation 
potential. Error bars indicate the quadratic sum of the error on the solar-system 
abundance and that on the cosmic-ray source abundance. The element symbols for 
the various points have been transferred directly below the point. 

Cosmic-ray source abundances in all four figures are derived principally from 
observations on HEAO-3, except for H and He which come from balloon observations 
(Webber, 1932; Webber and Lesniak, 1974). Two H points are plotted; H(R) uses H 
abundances at the same rigidity as the other elements; H(E), at the same energy per 
nucleon. For 6 < Z < 27, and for Z — 29 and 31 (Cu and Ca), the results are from the 
Danish-French experiment on HEAO-3 (Lund, 1984). For Z — 28 and 30 (Ni and Zn) 
and for Z > 32 the results are from the Heavy Nuclei Experiment (Israel, et a L, 1983). 
In each case experimental results were propagated back to the source in a standard 
leaky-box model. 

In figures 1 and 3 the solar system abundances are from Anders and Fbihara 
(1982). ^hese abundances are mainly from type Cl meteorites; except H, C, N, and O 
are from photospheric measurements, He is from the solar wind H/He ratio, Ne is from 
the solar wind Ne/Ar ratio and from astronomical measurements of extra-solar -system 
nebulae, and Ar, Kr, and Xe are interpolated from nearby elements. In figures 2 and 4 
the Cl meteorite abundances are replaced by C2 meteorite abundances (Mason, 1979) 
in forming the solar-system values. 

Figures 3 and 4 have the same ordinates as figures 1 and 2 respectively, but the 
abscissa it A/Q^, where A is the atomic weight and Qj-iq is the charge state the 
eleme.it would have after removal of all electrons whose ionisation potential is less 
tba.i 120eV. Thus the abscissa is an estimate of the mass-to-charge ratio which the 
element would have in a million-degree interstellar gas. (The value 120 eV corresponds 
to a temperature of 1.4 x 10° K, but the correlation displayed here is insensitive to the 
precise value of ionization potential selected.) 

The implications of these figures will be discussed at the conference. 
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INTERACTIONS OF HEAVY NUCLEI, Kr, Xe AND Ho, IN LIGHT TARGETS 
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C.J. Waddington, 1 W.R. Binns, 2 T.L. Garrard, 3 M.h'. Israel 2 
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Pasadena, CA 91125 

1. Introduction. Over the past few years, we have been analyzing the 

IIEA0-3 measurements of the abundances of ultra-heavy cosmic ray nuclei 
(Z>26) at earth. 1 In order to interpret these abundances in terms of a 
source composition, allowance must be made for the propagation of the 

nuclei in the interstellar medium. Vital to any calculation of the 

propagation is a knowledge of the total and partial interaction cross 
sections for these heavy nuclei on hydrogen. Until recently, data on 
such reactions have been scarce, and we have relied on the semi- 
empirical formalism of Silberberg and Tsao 2 to predict the partial cross 
sections. However, now that relativistic heavy ion beams are available 
at the LBL Bevalac, some of the cross sections of interest can be 
measured at energies close to those of the cosmic ray nuclei being 
observed . 

During a recent calibration at the Bevalac of an array similar to 
the HEA0-C3 UH-nuclei detector, we exposed targets of graphite (C), 
polyethylene (CH n ), and aluminum to five heavy ion beams ranging in 
charge (Z) from* 36 to 92. Total and partial charge changing cross 
sections for the various beam nuclei on hydrogen can be determined from 

the measured cross sections on C and C^, and will be applied to the 

propagation problem. The cross sections on A1 can be used to correct 

the abundances of UH cosmic 
rays observed in the HEAO 
C-3 detector for 

interactions in the 
detector itself. Table 1 
shows the combinations of 
beams and targets, as well 
as the number of events 
incident on the target for 
each run. The energies of 
each beam are also shown. 

Our preliminary results 
show that we achieved a 
charge resolution on the fragments that ranged from 0.21 charge units 
for Kr on A1 to 0.28 c.u. for Au on C, permitting unambiguous resolution 
of individual fragments. In this paper we report on the total cross 
sections for Kr on A1 , and total and partial cross sections for Xe and 
Ho on C, CH 2 and H. 

2. Experimental Setup. The detector consisted of an array of two front 

ion chambers, a target space, two rear ion chamber. 0 , followed by a Pilot 
425 Cherenkov counter. From signals in the front ion chambers, we find 


Table 

1 . 

Number 

of Events (x 

10 ) 

Energy 

Ki 

Xe 

Ho 

Au 

U 

(GeV/n): 1.5 

1.2 

1.1 

1.0 

0.9 

Targe'. 






C 



210 

130 

260 

60 

CH 2 

— 

330 

200 

VOO 

90 

A1 

90 

160 

190 

200 

— 

"Blank" 

35 

110 

120 

260 

40 


r: 


chat approximately 10Z of the nuclei incident on the detector do not 
have the nomii.al charge of the beam, and we eliminate these events from 
further analysis. Fragments produced in the target as well as beam 
nuclei surviving through the target are measured in the rear ion 
chambers and the Cherenkov counter. Scatter plots of the signals in 
these detectors show well resolved peaks for individual fragments. 

3. Total Cross Sections. The interaction mean free path can be found by 
counting the number of beam nuclei which survive through the target and 
our detector, and correcting for interactions in the detector itself. 
This correction is found from a "blank" or no target run. The total 
interaction cross section per nucleus is related to the mean free path 
by the following expression: 

0 (mb) = A^,/(6.02 x 10 4 ) A (g/cm^) 

where is the mean mass number of the target. Our results for the 
total cross sections are given in Table 2, along with values calculated 
using the formula from Westfall et al. 3 for charge changing cross 
sections, 0^. Although this formula was derived from data for nuclei 

with Zs26, it gives 
Table 2. Total Cross-Sections values which only 


3eam 

Target 

K g/ cm 2 ) 

<^mb) 

1 

n 

slowly deviate from 
those measured as Z 
increases, with o/o u = 
0.85 for Ho on H. 

Kr 

Al 

19.3 

2300±100 

2460 

.95 

Xe 

Xe 

C 

CH 2 

8.9 

5.1 

2240±80 

1510±50 

2460 

1670 

.91 

.91 

4. Partial Cross 

Sections. The numbers 

Xe 

H 

1.4 

1150±90 

1270 

.90 

of fragments produced 
in the target are 

Ho 

C 

7.8 

2560±70 

2760 

.93 

measured in the rear 

Ho 

CH- 

4.6 

1690±40 

1910 

.89 

ion chambers and the 

Ho 

L 

H 

1.3 

1260±75 

1490 

.85 

Cherenkov counter. 

A histogram of the 


Cherenkov signal for events consistent with their being fragments is 
shown in Fig. 1. The charge resolution for this particular run is 0.23 
charge units. 



22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 


CHARGE CHANGE, AZ 

Fig. 1. Cherenkov histogram for Xe nuclei on Carbon. 

Several corrections must be applied to the yields obtained from 
these histograms. First, there is a correction for events observed in 
the blank run which are also consistent with their being fragments. 
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These events are due to interactions of the beam in the matter between 
the Cherenkov counter and the ion chambers. Also, for a given Z, there 
is a background due to the fragments of charge Z + 1 making a AZ = 1 
interaction in the Cherenkov counter. Both of these corrections are 
small, being less than 2%. In addition, there must be a correction for 
absorption in the detector. We have done this by applying an 
exponential absorption law, using the Westfall et al. 3 charge changing 
cross section, scaled to our measurements, to calculate the necessary 
mean free paths. The numbers resulting from this correction are the 
numbers of each fragment exiting the target. To obtain the partial 
cross sections from these numbers we need to correct for multiple 
interactions in the target. The targets used were approximately 0.25 of 
a mean free path, and were chosen as a compromise between being thick 
enough to produce a reasonable number of interactions, yet thin enough 
so as not to degrade the charge resolution of the Cherenkov due to Lr.e 
energy spread of the fragments. We have used a slab propagation program 
to do this thick target correction. 

Table 3 lists the partial cross sections of Xenon and Holmium on C, 
CH^, and H. The hydrogen cross sections are derived from the C and CH^ 
cross sections per nucleus by a subtraction procedure: 

° H = 1/2 (3 °ch 2 ~ °C ) 

Also given in Table 3 are the values predicted 4 for Xe and Ho on H. 
Fig. 2 shows the ratios of our values and those predicted, as a function 
of AZ for Xe and Ho. The errors shown are the statistical errors on the 
target and blank runs, combined with the errors due to the top of 
detector correction. Also shown is a fit to previously measured ratios 
reported for Au nuclei, 5 showing distinctively different behavior. 

Table 3 


a 

Xe on C 

Xe on CH 2 

Xe on H 

S & T 

Ho on C 

Ho on CH^ 

Ho on H 

S & T 

i 

249±21 

230*9 

220*17 

257 

343*35 

270*26 

234*42 

267 

2 

128*5 

137*4 

141*6 

169 

138*5 

133*4 

131*6 

257 

3 

105*4 

104*3 

104*5 

106 

103*4 

117*3 

124*5 

166 

4 

72*4 

94*3 

105*5 

106 

81*4 

97*3 

105*5 

114 

5 

73*4 

79*3 

82*5 

61 

73*3 

89*3 

97*4 

95 

6 

63*3 

73*2 

77*4 

64 

67*3 

78*2 

84*4 

81 

7 

53*3 

65*2 

71*4 

50 

60*3 

67*2 

71*4 

66 

8 

53*3 

53*2 

53*1 

47 

53*3 

57*2 

59*3 

71 

9 

45*3 

47*2 

48*3 

33 

45*3 

51*2 

54*3 

53 

10 

43*3 

43*2 

43*3 

31 

46*3 

49*2 

51*3 

51 

1 1 

40*2 

37*2 

35*3 

21 

43*3 

40*2 

38*3 

34 

12 

35*2 

33*2 

32*3 

18 

34*2 

32*7 

31*3 

32 

13 

37*2 

31*2 

28*3 

14 

38*2 

29*2 

24*2 

18 

14 

37*2 

27*2 

23*3 

14 

33*2 

25*1 

20*2 

14 

15 

33*2 

20*1 

14*2 

9 

33*2 

22*1 

17*2 

12 

16 

33±2 

19*1 

12*2 

8 

32*2 

19*1 

13*2 

8 

17 

30*2 

17*1 

1012 

6 

34*2 

13*1 

10*2 

7 

18 

32*2 

16*1 

7*2 

6 

32*2 

15*1 

7*2 

5 

19 

31±2 

15*1 

6*2 

5 

27*2 

13*1 

6*2 

5 

20 

30*2 

13*1 

5*2 

6 

25*2 

15*1 

10*2 

5 

21 

27*2 

13*1 

6*2 

5 

— 

— 

— 

— 
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Fig. 2. Ratio of observed to predicted cross-sections as a function of 
AZ. 


Examination of these data and of the earlier Au-data from Brewster 
et al. (1983) shows that they can be well represented by universal 
curves if expressed as dO/O^, versus AZ. We find that for heavy targets, 
carbon and aluminum, do/O- = a (AZ) where a and m are closely similar 
constants for all non-hydrogenous targets and projectiles. Similarly 
for a hydrogen target, do/o T = b exp (-nAZ), where b and n are closely 
similar constants for all studied projectiles. The polyethylene targets 
also show a similar exponential dependence. 
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ABSTRACT 

As part of a recent calibration at the LBL Bevalac for the Heavy 
Nuclei Experiment on HEAO-3, we have compared the response of a set 
of laboratory ionization chambers to beams of 26 Fe, ^Kr, j 4 Xe, 67 Ho, and 
79 AU nuclei at maximum energies ranging from 1666 MeV/amu for Fe to 
1049 MeV/amu for Au. The response of these chambers shows a 
significant deviation from the expected energy dependence, but only a 
slight deviation from Z 2 scaling. 


1. Introduction 

Gas filled ionization chambers were used on the Heavy Nuclei Experiment (HNE) 
on HEAO-3 (Binns et al., 1981). The response of such chambers is expected to be pro- 
portional to the energy deposited by the particle traversing them. At low energies this 
energy deposit is simply the ionization energy loss, while at high energies energetic 
knockon electrons are able to escape from the chamber, reducing the energy deposit. 

To first order the ionization energy loss scales as the square of the particle charge Z, 
hov’ever at high Z this assumption breaks down. A more complete expression is given by 
Ahlen (1980, 1982), and predicts an energy loss rising slightly faster than Z 2 . Such effects 
are important when identifying ultraheavy elements. 

We have performed two calibrations of ion chambers at the LBL Bevalac using 
beams ranging finm 23 Mn to 79 AU. Th' first, in 1982, was done with a prototype of the 
HNE ion chamber module which was essentially identical to that used in flight. Thus 
those data, reported in Garrard et al., 1983, are directly applicable ro our flight experience 
at the energies calibrated. The second calibration, in 1984, used lab chambers which were 
made of thinner and more uniform materials, permitting better resolution and better 
knowledge of the beam energy in each ion chamber, at the cost of less direct relevance to 
the flight data. Figure 1 is a schematic drawing of the 1984 detector. 

Particles entering the 1984 detector traversed ~-0.1 g cm -2 of mylar in the upstream 
window, rather than the — 1 g cm * 2 of aluminum honeycomb in the flight prototype; thus 
the energy loss in the window is much smaller and more uniform. Also, in the 1984 cali- 
bration the beam energy was measured with a magnetic spectrometer after being degraded 
to the calibration energy, rather than being calculated from an energy loss model. 

The lab ion chambers had aluminized mylar electrodes (0.8 mg cm* 2 ) rather than 
aluminum screenwire (10 mil diameter, 62.5 rail spacing); thus the production and 
absorption of knockons is much more uniform. A Monte Carlo model of knockon 
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production correctly predicts the degradation in resolution caused by non-uniform pro- 
duction of knockons in the screen wire electrodes. This resolution degradation in the 
flight chambers tends to mask the rela ively subtle deviations from Z 2 scaling. 



Figure 1. Schematic drawing of the 1984 detector, showing the six ion chambers. 

In 1982 the ultraheavy capabilities of the Bevalac were new and we calibrated only 
on beams of — 1700 MeV/amu 23 Mn and —1000 MeV/amu 79AU. The 1984 calibration 
used beams of 26 Fc, *Kr, ^Xe, 67 Ho, and t*Au at maximum energies ranging from 1666 
MeV/cmu for Fe to 1049 MeV/amu for Au. 

2. Results of the 1984 Calibration 

Figure 2 shows the response of chambers 1, 5, and 6 to 26 Fe nuclei as a function of 
the energy at the midplane of the appropriate chamber, and compares their signals to the 
calculated dE/dx, arbitrarily normalized at 500 MeV/amu (requiring 27.9 eV per ion pair 
in the P-10 gas used (90% argon, 10% methane)). It is apparent that the signals fall below 
that predicted by dE/dx at energies above 700 MeV/amu. This loss of signal is somewhat 
surprising since at these energies we would .xpect knockons escaping from the exit win- 
dow to be in equilibrium with those arriving from above, particularly for chambers 5 and 
6 which have —2 g cm -2 of upstream material. However, some of the decrease in 
observed signal may be due to knockons escaping from the .ides of the chambers. 

By interpolating to a particular energy we can construct a plot of signal versus Z at 
that energy. At low energies, the heaviest nuclei have an effective charge, 
Zrf- Zll— exp(-130pz _2/3 )], due to electron capture (Pierce and Blann, 1968). Figure 3 
shows the pulse heights, scaled down by Z^ ff . at four energies for Z - 26-79, using ion 
chambers 1, 5, and 6. The uranium data have not been included because the charge state 
in the magnetic spectrometer was uncertain for those beams whose energy had been 
degraded significantly. The straight lines represent a linear fit to the data, and it is 
apparent that there is a small negative non-Z 2 effect. The charge of an 82 Pb nucleus would 
be underestimated by about 0.5 charge units at these energies, in contrast with the charge 
overestimate of +3 charge units observed in the calibration of the flight chambers. 
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3. Conclusions 

Although the non-Z 2 effects in these chambers differ from those obse.ved in the pro- 
totype flight chambers, the assumption of Z 2 scaling is still not seriously in error. We also 
.iote that our published abundances above charge 50 have used primarily the Cerenkov 
selector to assign charges, and are unaffected by small non-Z 2 effects in the ion chambers. 


Since the two calibrations differ, the ionization response to energy loss must be sen- 
sitive to details of the mass distribution above, below, and within the chambers. As a 
result we have used the flight data to directly determine both the energy dependence and 
effective non-Z 2 correction (Jones et al., (1985, OG 4.1-8) °.nd Newport et al., (1985, OG 
4.4-5)). 
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